Discrete vortex solitons and VT symmetry 
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We study the effect of lifting the degeneracy of vortex modes with a PT symmetric defect, using discrete 
vortices in a circular array of nonlinear waveguides as an example. When the defect is introduced, the 
degenerate linear vortex modes spontaneously break PT symmetry and acquire complex eigenvalues, but 
nonlinear propagating modes with real propagation constants can still exist. The stability of nonlinear modes 
depends on both the magnitude and the sign of the vortex charge, thus PT symmetric systems offer new 
mechanisms to control discrete vortices. © 2013 Optical Society of America 
OCIS codes: 230.7370,080.6755,190.3270 
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Nonlinear periodic media allow for efficient localiza- 
tion and routing of light signals in the form of discrete 
optical solitons [1]. Alongside usual signals encoded in a 
soliton's amplitude, the phase of a vortex soliton also car- 
ries a quantized bit of information, its topological charge 
m (TC) [2]. The absolute value of TC, |m|, is the phase 
winding number, in units of 27r, around the vortex ori- 
gin, and the sign of TC determines the direction of power 
flow, or vortex handedness. Periodic media support sta- 
ble discrete vortex solitons [3] and allow for robust con- 
trol over vortex TC by switching between integer values 
m o — m, known as charge flipping [4-6]. 

In systems with discrete rotational symmetry, such as 
a ring of N identical waveguides [6], vortex and anti- 
vortex modes with phase ~ eyi'p{i2T:mn/N) are degen- 
erate, here n = 1, 2 ... A" is the waveguide number. The 
charge-flipping transformation m o —m can be seen 
as a complex conjugation, or time-reversal operation. It 
is interesting to explore the additional possibilities for 
TC control which can be offered by dissipative parity- 
time {'PI') symmetric systems [7], which have pure real 
spectra below the T'T'-symmetry breaking point. Recent 
studies of infinite chains and discrete rings of coupled op- 
tical waveguides with gain and loss have demonstrated 
their direct relevance to applications [8-12]. The exis- 
tence of propagating modes with real spectra in such 
arrays readily suggests the possibility of the existence of 
propagating vortices. Broken T symmetry implies that 
vortices with opposite charges can behave differently, i.e. 
their degeneracy is lifted. 

In this Letter we show that nonlinear propagating vor- 
tex modes can still exist even when the VT symmetry 
breaking threshold is reached for the linear system. We 
explore interesting consequences of lifted vortex degener- 
acy. The onset of modulational instability becomes sen- 
sitive to both the magnitude and sign of vortex charge, 
thus expanding the available toolbox for TC control. The 
sensitivity to the sign of vortex charge is impossible in 
systems respecting 7" symmetry. 



We consider a circular array of N waveguides with 
Kerr nonlinearity, shown in Fig. 1(a), with gain and 
loss located at waveguides 1 and N, coupling C between 
them, and nearest neighbor coupling normalized to unity 
along the rest of the ring (the difference in the coupling 
constants can be introduced by variation of the distances 
between the respective waveguides). 
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Here En{z) {n = 1...N) is a dimensionless electric field 
in the n— th waveguide, En+N = E„^ and S = ±1 is 
the nonlinear coefficient. Stationary modes take the form 
En{z) — >■ ii'n exp(i/3z), where (3 is the propagation con- 
stant. Aspects of the linear, large N limit of Eq. (1) were 
previously considered in Refs. [12,13]. In discrete conser- 
vative systems the TC is m = ^ ^^^^ Arg[i?*i?„+i], 
and measures the phase winding along the contour. 
Jn = 2Im(ii^*i?„+i) is the energy fiow from site n to 
site n + 1 [6] . The usual assumption is that nonzero TC 
(m > 0) m < corresponds to the (anti)clockwise cir- 
culation of phase and energy around the contour. This 
definition breaks down for even N at the band edge 
|to| = A/2, which are multipole modes En ~ (— 1)" 
without energy circulation (ie. not a vortex). 

The equation for stationary solutions in ring-type sys- 
tems can be cast in the general form 

-I3E + HE + 5F{\E\^)E^Q, H = Ha + i"fHi (2) 

where E = {Ei, ...^En)^ is a column vector {T stands 
for transposition) and the nonlinearity is given by the 
diagonal matrix F{\E\'^) =diag(|i?i|^, [iJjvl^)- Linear 
operator Hq is a matrix describing the array without 
dissipation and Hi describes the losses. One ensures that 
the commutator [Hq, T'] = and HiP = —VHi with the 
"parity" inversion matrix V having only anti-diagonal 



nonzero elements, Vi 



and that [i/o,i, T] 
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Fig. 1. (Color online) (a) Schematic of iV-site ring with 
gain (+) and loss (-) at waveguides 1 and N. Phase cir- 
culation direction of vortices with m > and m < 
is indicated by anti-clockwise (red) and clockwise (blue) 
arrows, (b) /? vs. C for a conservative (7 = 0) = 4 ring. 
Degenerate vortex modes only occur at the intersection 
marked by the black circle, (c) and (d) Linear spectrum 
for = 3 and = 4 rings with 7 ~ 0.2. Modes with 
TCs -1-1 (-1) are shown in red (blue). 



0, where T is the operator of complex conjugation. Thus 
[H, VT] = 0, i.e. H is PT-symmetric. 

First, we show that if Hq supports at least one de- 
generate pair of vortex modes E± belonging to a double 
degenerate eigenvalue /3o, i-C. HqE± ~ (3oE±, and Hi is 
symmetric, such as in Eq. (1), the VT symmetry break- 
ing threshold 'jth, i-C. the value of 7 above which the 
spectrum has complex eigenvalues, is zero: jth = 0. In- 
deed, from the above properties of Hq it follows that 
one can choose E± = VEzfi = TE^ = E^. Hence 
EW =E+ + E. and E''^) =i(E+-E^) are the two real 
eigenstates of Hq and VE^^'' = {-iy+^E'^^'> {j = 1,2). 
Then, from the Theorem 2.1 of Ref. [14] follows that H 
has a pair of complex eigenvalues for arbitrarily small 7. 

The symmetry breaking threshold of Eq. (1) at C ^ 1 
is non-zero [12]. Hence, there are no vortex eigenstates 
of Hq. Instead, vortex eigenstates appear at nonzero 
7 < jth- Moreover, the VT operator does not change 
TC, i.e. if E is an eigenmode of H with the charge m, 
HE = (3E at 7 < jth, then VTE is also an eigenmode 
with the same /3 and m, and thus one can set VTE = E. 
In other words one can consider linear PT-symmetric so- 
lutions. In Fig. 1(b) we show the spectrum for an A'^ = 4 
ring when 7 = 0. Vortex modes only exist at C = 1, and 
they are degenerate. The effect of nonzero 7 is shown in 
Figs. l(c,d) for odd A^ = 3 with 7t,, = {C^ -3C^/^ + 2)^/^ 
and even A^ = 4 with jth = \C — 1|. The degeneracy be- 
comes a pair of exceptional points, and branches with 
nonzero topological charge |m| = 1 appear. The break- 
ing of T symmetry manifests itself through the separa- 
tion of different charges into distinct branches. Curiously, 
nonzero m here does not imply nonzero vorticity: the 




Fig. 2. (Color online) Stable (unstable) nonlinear modes 
of A^ = 3 ring are shown with solid (dashed) lines. TCs 
are indicated next to the curves, m = +1 in red (purple), 
m = — 1 in blue (brown), and m = in black (grey). 
Parameters are: (a) C = 0.6,7 = 0-2 < 1th', (b) C = 
0.6,7 = 0.6 > jth; (c) C = 1.3,7 = 0.2 < jth] and (d) 
C = 1.3,7 = 0.6 > 7*;,. 



multipole (lowest (3) mode in Fig. 1(d) has jti = 1, but 
it is not a vortex: there is phase winding, but its phase 
does not increase monotonically, resulting in flow from 
the gain site to the dissipative site along both paths. In 
contrast the |m| = 1 modes which vanish as C — >■ 1 are 
true vortices with energy circulation (the sign of J„ is 
the same for all n). 

Since the nonlinearity in Eq. (2) has the symmetry 

VT (^F{\E\'^)E^ = F{\E\'^)E there exist [11] nonlinear 
propagating vortex modes bifurcating from each of the 
linear 7^7~-symmetric vortices. These modes do not ex- 
haust all possible soliton solutions with real propagation 
constants which can exist even when VT symmetry is 
broken in the linear regime [11]. In Fig. 2 we present 
families P{I3) of nonlinear modes obtained using New- 
ton's method (here P = J2n=i the total power). 
Note that for 7 ^ the total power is not conserved 
except for stationary modes. We consider below the sim- 
plest A^ = 3 ring in detail, then discuss properties of 
larger A^ systems. 

Fig. 2(a) shows the typical A^ = 3 soliton family for 
C < 1 with 7 < 7t/i. We observe the bifurcation of 
m = 1 vortices from linear modes at P = discussed 
in Fig. 1(c). As the power is increased, the energy flow 
between Ei and En in the lower branch decreases and 
eventually changes sign, destroying the vortex, and the 
charge becomes zero. No to = — 1 modes arc found in the 
parameter range scanned: there are no such linear modes, 
and the broken T symmetry suppresses their saddle-node 
bifurcation. 

When 7 is increased beyond jth in Fig. 2(b), the lin- 
ear modes merge and annihilate, and pairs of m = -1-1 
vortices now appear at a saddle-node bifurcation (notice 
that the unstable branch is not visible on the scale of the 
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figure). One branch quickly loses its charge and vortic- 
ity; this happens at lower powers as 7 is increased. Once 
again, no m = — 1 modes are found. 

The situation for C > 1, [Figs. 2(c,d)] is different. 
Now nonlinear m ~ —1 modes exist below ^th, because 
they can bifurcate from linear modes. In this case one 
is stable at high power, while the other is unstable, but 
both maintain their vorticity. A pair of m = +1 vortices 
is still created at a saddle node bifurcation, with behavior 
similar to Figs. 2(a,b); one branch loses its charge and 
vorticity, while the other remains stable. Above the VT 
breaking threshold, the m = — 1 modes are destroyed, 
while the m = +1 saddle node bifurcation remains. 

Evidently, the VT symmetric defect in combination 
with the asymmetric coupling C determines the charge of 
linear vortex modes - the energy flows from the site with 
gain to the site with loss along the path with stronger 
coupling, and back via the path with weaker coupling. 
In contrast, the charge of nonlinear modes above 74/1 
is mainly sensitive to the distribution of gain and loss. 
These charges are opposite when C > 1. 

As N is increased modes with larger \m\ < N/2 ap- 
pear. In conservative systems with focusing nonlinearity 
the modes with higher \m\ are stable, while lower charged 
vortices suffer instabilities above a critical power [6] . This 
behavior holds in VT symmetric rings, but the sign of 
the charge also plays an important role. For even N there 
is an additional nonlinear branch bifurcating from the 
linear multipole mode, which formally has the highest 
TC \m\ = N/2 — 1 supported by the ring, but it does not 
form a vortex. 

As an example of nonlinear vortex dynamics, we con- 
sider a hexagonal (A^ = 6) coupler, with C = 1 and 
6 = 1. When 7 = 0, the symmetric vortex modes of 
charge m take the form £'„ = Aexp{2mmn/N + i(3z) [6] 
and their stability is independent of the vorticity direc- 
tion, sign(m) . The situation is different with T symme- 
try broken, as we demonstrate in Fig. 3 for 7 = 0.2 and 
A = 0.5, by contrasting the propagation of m = ±2 
vortices. The m = +2 perturbed vortex oscillates sta- 
bly about a nonlinear vortex mode, and its charge is 
conserved. In contrast, the m = —2 input experiences 
VT breaking and its power grows exponentially. These 
dynamics are further illustrated by obtaining vortex 
lines from the array supermode (cf. Ref. [6]), shown in 
Figs. 3(c,d). We see the stable oscillation of a pair of 
charge 1 vortices for the m = +2 input, while the op- 
posite input charge m = — 2 results in highly irregular 
vortex dynamics. Inputs with \m\ = 1 arc also unsta- 
ble. With repulsive nonlinearity 6 = —1 the stability 
inverted: the m = 1 vortex is stable, while all others are 
unstable. In this example the nonlinearity is essential: 
VT breaking occurs in the linear limit, and vortex in- 
puts of any charge are decomposed into amplified and 
attenuated components. 

In conclusion, we have studied discrete vortices in A^- 
site rings of coupled nonlinear waveguides with a VT 
symmetric defect. We have shown that by breaking T 
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Fig. 3. (Color online) Charge selectivity in an array with 
iV 6, C = 1, 7 = 0.2, and S = 1. The total power 
(a) and TC (b) are shown vs. propagation distance z for 
m = -1-2 (black, dashed) and m = — 2 (gray, solid) inputs. 
Vortex (red) and antivortex (blue) lines are shown in (c) 
for input m = -\-2 and in (d) for m = —2. 



symmetry, the existence, stability and dynamics of non- 
linear vortex modes become sensitive to the sign of their 
charge, offering an additional degree of freedom for all- 
optical control of discrete vortices. 
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